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ABSTRACT 

We report mid-IR interferometric measurements with ~ 



10 mas resolution, which resolve the warm 



(T = 285±§g K) thermal emission at the center of NGC 4151. Using pairs of VLT 8.2 m telescopes 
with MIDI and by comparing the data to a Gaussian model, we determined the diameter of the dust 
emission region, albeit only along one position angle, to be (2.0 ± 0.4) pc (FWHM). This is the first 
size and temperature estimate for the nuclear warm dust distribution in a Seyfert 1 galaxy. The 
parameters found are comparable to those in Seyfert 2 galaxies, thus providing direct support for 
the unified model. Using simple analytic temperature distributions, we find that the mid- infrared 
emission is probably not the smooth continuation of the hot nuclear source that is marginally resolved 
with K band intcrferometry. We also detected weak excess emission around 10.5 /*m in our shorter 
baseline observation, possibly indicating that silicate emission is extended to the parsec scale. 
Subject headings: techniques: high angular resolution — galaxies: active — galaxies: Seyfert 



1. INTRODUCTION 
In the standard model of Active Galactic 



(AGNs) ( |Antonucci||1993| | Urry fc Padovani 
central engine (Black Hole 



1995 



Nuclei 
, the 
road 



accretion dis 

Line Region) is embedded in a dusty torus. Galaxies 
where we see the center only through this torus are called 
"type 2" whereas ones where we receive direct radiation 
from the Broad Line Region are referred to as "type 1" . 
For a long time, there was only indirect evidence for the 
existence of the enshrouding dust, until VLTI observa- 
tions with MIDI actually resolved the parsec-scale AGN 
heated dust structures in two Seyfert 2 galaxies (Sy 2) . 



( |Jaffe et al.|2004||Tristram et al. 12007] |Raban et al.|2009| ) 

Further observations proved the existence of parsec 
siz ed dust structures in several Seyfert 1 galaxies (Sy 
1) ( Tristram et al.|2009 l, but so far, no Sy Is have been 
observed in sufficient detail to test the central premise 
of the unified models: type 1 and type 2 dust distribu- 
tions are identical and the observed differences are only 
due to differences in orientation with respect to the line 
of sight. These predictions can now be tested with the 
MIDI observations of NGC 4151 reported here. 

At a distance D = (14 ± 1) Mpc (i.e. 1 mas ~ 
0.068 pc) 1 , NGC 4151 is the closest and brightest type 
1 galaxy (classification: Seyfert 1.5). It is also one of 
the most variable Seyfert galaxies: Th e UV conti nuum 
flux varies on scales of days and weeks ([Ulrich 2000J and 
the reverberation time to the hot dust on the sub-pars ec 
scale varies on yearly timescales (Koshida et al. 2009). 

Sy 1 galaxi es have been observe d previously with MIDI 
(NGC 378 3, |Beckert et al]|2008| , (NGC 7469, |Tristram 
et al.|2009 1 , but 1NGC 4151 is the first multi-baseline case 



*BASED ON OBSERVATIONS COLLECTED AT THE EURO- 
PEAN SOUTHERN OBSERVATORY, CHILE, PROGRAMME 
NUMBER 081.B-0092(A). 

from the NASA Extragalactic Database: 

|http:/ / nedwww. ipac.caltech.edu distance from redshift with 
Ho=73 km/s/Mpc; other estimates range to 20 Mpc 



where the size of the nuclear dust distribution is clearly 
indicated. 

al. (2003) reported near-IR, 2.2/im, observa- 
the Keck interferometer . They find that the 



Swain et 
tions with the 



majority of the K band emission comes from a largely 
unresolved source of < 0.1 pc in diameter. Based on this 
small size, they argued that the K band emission arises 
in the central hot accretion disk. They note, however, 
that their result is also consistent with very hot dust at 
the sublimation radius. This view is s upported by the 
K ban d reverberation measurements of Minczaki et al.l 
(2004) who find a l ag time corresponding to ~ 0.04 pc. 



Kitfel et al. ( 2009 1 modelled the near-infrared spectrum 
of this source and found it to be composed of a powerlaw 
accretion-disk spectrum and a component likely arising 
from hot dust (T = 1285 K). The hot dust component 
dominates at A > 1.3/im, consistent with the interpreta- 
tion of hot dust emitting at 2/im. They measured a K 

band flux of ~ 65 mJy. 

In the mid-infrared, Radomski et al. (2003) presented 
images of NGC 4151 at I0.8/zm with an aperture of 4.5" 
and a spatial resolution of ~550 mas (35 pc). They find 
that the majority (73%) of the N band flux comes from 
an unresolved point source with a size < 35pc, and the 
rest is extended emission from the narrow line region. 

In this letter we report new mid-infrared interferomet- 
ric observations of NGC 4151 which clearly resolve a ther- 
mal structure. 

2. INSTRUMENT, OBSERVATIONS AND DATA 
REDUCTION 

Observations were performed in the N band (A rs 
13/im) with t he MID-infrare d Inte rferometric in- 



strument (MIDI, |Leinert et al.| |2003[) at the ESO 
Very Large Telescope Interferometer (VLTI) on Cerro 
Paranal, Chile, using pairs of 8.2 m Unit Telescopes 
(UTs). The observables with MIDI are the single-dish 
spectra and the correlated flux spectrum that is ob- 



2 



tained from the interference pattern generated by the two 
beams. The spectra were dispersed with a NaCl prism 
with R = X/5X ~ 30. 

Observations were taken in the nights of April 21 and 
23, 2008 with projected baseline length BL of 61 m and 
89 m at position angles of 103° and 81° respectively. 
These provide effective spatial resolutions 2 (X/3BL) at 
10.3/mi of 11 and 7 milliarcseconds (mas), respectively. 
The calibrators HD 133582 and HD 94336 were selected 
to be very close in airmass with A(secz) < 0.15. This 

is especially important for NGC 4151 (DEC 1-40°) 

which, at Paranal, never rises higher than ~ 25° (sec z > 
2.3) above the horizon. The northern declination of the 
source also essentially limits the projected baselines and 
fringe patterns to East- West orientation at Paranal. The 
N band spectrum of HD 133582 (K2III) was taken to fol- 
low a Rayleigh- Jeans law, wh i le tha t of HD 94336 (Mill) 



was taken from Cohen et al. ( |1999| . 



Data reduction was performed with the interferometric 
data reduction software MIDI Inter active Ana lysis and 
Expert Work Station (MIA+EWS, [Ja^|[2004| , specifi- 
cally adapted for the analysis of low S/iN sources. 

3. RESULTS AND MODELLING 
3.1. Single-dish spectrum 

The resulting single-dish spectrum with an effective 
aperture of ~ 300 mas is shown in Fig. [I] together with 
a Spitzer IRS spectru m for this source (obs erved 8 Apr 
2004; aperture ~ 3.5", |Weedman et al.|2005l . The higher 
Spitzer flux most likely indicates emission from the Nar- 
row Line Region (cf. |Radomski et aT||2003[ ). 

The color temperature of the spectrum is (2851 5 q) K. 
Additionally, we detect the [Ne II] 12.81 /im line com- 
monly seen in star forming regions. The [S IV] (10.51 
/mi) line, clearly seen in the Spitzer spectrum, is not sig- 
nificant in the MIDI spectrum. The errors in the MIDI 
spectrum arise from incomplete thermal background sub- 
traction and hence rise steeply with increasing wave- 
length. 

3.2. Correlated spectra 

The two correlated flux spectra are shown in Fig. [2j 
The correlated flux observed with BL = 61m rises 
from (0.13 ± 0.01) Jy at 8.5 /mi to (0.43 ± 0.05) Jy at 
12.5 /mi. The second spectrum (BL = 89m), rises from 
(0.12 ± 0.02) Jy at 8.5 /mi to (0.30 ± 0.05) Jy at 12.5 
/im. Correlated flux uncertainties arise primarily from 
background photon noise and increase with wavelength 
but are smaller than the single dish errors. 

The correlated flux on the shorter baseline (the one 
that has a higher flux) shows a broad "bump" between 
9 /mi and 12 /im th at w e interpret as a silicate emission 
feature (see section [33| . No [Ne II] emission is seen, in- 
dicating that this arises on a scale that is fully resolved 
by the interferometer (i.e. > 20 mas ~ 1.3 pc). The fact 
that the correlated flux is lower on the longer baseline is 
a clear sign that the source is resolved by the interferom- 
eter. 

2 Applying Raylcigh's criterion to an interferometer leads to 
a resolution of X/2BL. However, given sufficient S/N, one can 
already distinguish models at lower resolutions, in our case at 
~ X/3BL. 
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Fig. 1. — Si ngle-dish spectra for NGC 4151. Spitzer spectrum 
(3.5", dotted, |Weedman et al.|2005| , MIDI spectrum (0.3", black 
line with error bars). The MIDI single-dish errors were taken as the 
error of the mean from five observations. Also plotted are black- 
body emission curves at T = 235,285,310 K (green, red, blue). 
The region of atmospheric ozone absorption, between 9.3 and 10 
fim (hatched), is uncertain and not taken into account for the later 
analysis. 
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Fig. 2. — Smoothed (AA ~ 0.3/im) correlated MIDI spectra at 
two different East-West projected baselines: 61m (red, average 
of two fringe track observations) and 89m (blue, average of three 
fringe track observations). The errors are the errors of the mean 
of the individual observations. The region of atmospheric ozone 
absorption, between 9.3 and 10 /im (hatched), is uncertain and 
not taken into account for the later analysis, as are the regions 
A < 8.3/xm and A > 12.7/im (hatched), which have very low S/N. 

Since the resolution of the interferometer 9 m i n ~ 
X/'iBL changes with wavelength, the correlated flux re- 
flects both the source spectrum and the source structure. 
It is not possible to draw any conclusions from the spec- 
tral slope of a correlated flux spectrum without assuming 
a source geometry. 

3.3. A possible silicate emission feature 

Although the silicate absorption feature has often been 
detected in type 2 nuclei, the emissio n feature, predicted 



for type 1 nuclei by torus models (e.g. Pier & Krolik 1992 



Schartmann et al.||2005 ), has not been detec ted except 
in a handful of objects, most of them quasars (|Hao et al. 



2005| |Weedman et al.||2005| |Buchanan et al.||20D5| r 
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The silicate emission feature as s een in v 2 F v . In the 
Spitzer spectrum (black, Wecdman ct al. 20 05]) , a broad weak emis- 
sion feature at around II /im is clearly seen. The three MIDI spec- 
tra (green: singlc-dish, red: 61m baseline, blue: 89m baseline) are 
plotted with offsets. The grey line shows a 285 K blackbody. The 
region between 9.3/xm and 10 /im is hard to calibrate in the MIDI 
spectra due to the atmospheric ozone feature and has been greyed 
to not mislead the eye. 



As noted by Weedman et al. (20051 and 



Buchanan 

et al.| ( [2006| , Spitzer spectra (with an aperture of ~ 3.5") 
of JNGC 4151 show weak excess emission at 10 /xm and 18 
^m that is most easily seen when plotted as v 2 F v with a 
sufficiently large wavelength range. In Fig. [3] we plotted 
our spectra (that are limited to the atmospheric N band) 
in such a way together with a Spitzer spectrum. 

The hypothesis of silicate emission is not inconsistent 
with our single dish spectrum (aperture ~ 0.3") but since 
this spectrum suffers from incomplete background sub- 
traction, especially at longer wavelengths, it is probably 
hidden in the resulting uncertainties. The emission fea- 
ture seems to be most prominent in our observations on 
the 61m baseline observation and is clearly not detected 
on the 89m baseline. 

This suggests that silicate emission is distributed over 
~ 1 pc as derived from the 10.5 /im Gaussian model (see 



3.4. Simple Gaussian model 

With the limited baselines available from Par anal, it is 
not possible to reconstruct an image of NGC 4151 from 
our data. Instead we consider simple model distributions 
of the emission on the sky and compare the predicted in- 
terferometric and single dish spectra with our measure- 
ments in order to fix parameters in such a model. We 
chose a model containing an unresolved point source (flux 
F p ) and an extended Gaussian distribution (flux normal- 
ization F g , FWHM 9). Although we might expect the 
mid-IR brightness distribution in Sy 1 galaxies to have a 
hole in the middle, neither observed nor modelled fluxes 
change as long as the hole, whose radius is determined by 
the sublimation radius of the dust, is unresolved. This is 
certainly the case in NGC 4151 (see sectionR}. An upper 
limit to the size of our point source is givenby the effec- 
tive resolution of the interferometer; this corresponds to 
a diameter of ~ X/'dBL ~ 7 mas, i.e. a radius of ~ 0.2 
pc, at 10.3 /urn. 




2x10" 4x10" 6x10° 8x10" 
Spatial Frequency [fringe cycles / red] 

Fig. 4. — Gaussian model for NGC 4151. On the upper x axis 
the resolution (i.e. sensitivity to model parameters) of the interfer- 
ometer X/3BL is given. The three curves correspond to uniquely 
determined models at 8.5 fim (blue), 10.3 £im (green) and 12.5 £im 
(red). Data points and errors are taken from the single-dish (i.e. 
"0-baseline" , Fig. [I| and correlated fluxes (the remaining two data 
points per wavelength, Fig. [2j- The many thin lines represent the 
Monte-Carlo-like error determination for the Gauss model. For 
details see text. 



TABLE 1 

Parameters for the Gaussian models (see text). 



A [/im] 


F g [Jy] 


9 [mas] (diameter [pc]) 


F P [Jy] 


8.5 ±0.2 
10.3 ±0.25 
12.5 ±0.3 


0.41 ± 0.10 
0.70 ±0.16 
1.03 ±0.30 


29+ a _ 6 (2.0+ u ^_o. 4 ) 
23 ±4 (1.5 ±0.3) 
32 ±6 (2.1 ±0.4) 


0.119 ±0.016 
0.194 ± 0.061 
0.290 ± 0.070 



Because of the East-West baseline orientation, the 
North-South distribution of emission is undetermined. 
Equivalcntly, we assume the source to be circularly sym- 
metric on the sky. 

The model correlated flux density at wavelength A is 
then given by F V (X) = F p (X) + F g (X) ■ exp(-(B£/A ■ 
it/2 ■ 6>) 2 /ln(2)). With the given data points this model 
is uniquely determined. We calculated the parameters of 
such a model separately at 8.5, 10.3 and 12.5 where 
we are safely away from the regions of very low signal to 
noise (at the edges of the N band) and the ozone feature. 
At 10.3 ^m the parameter values may be affected by the 
silicate feature. The modelled visibilities are shown in 
Fig. [4] and the resulting parameters are given in Table [l] 

The errors of these parameters were estimated from 
a Monte-Carlo simulation in which we randomly placed 
10000 measurements in a Gaussian distribution around 
the measured value with the a as determined from our 
data reduction. The parameter errors are then given by 
the standard deviations of the resulting model solutions 
to the simulated data. 

4. DISCUSSION 

4.1. The extended source and the Sy 1 / Sy 2 paradigm 

In the strictest version of unified models, we expect for 
both a Sy 1 and a Sy 2 galaxy an extended dust struc- 
ture with the same size, morphology and temperature 



4 



distribution (at a given UV luminosity Ltjv)- In less 



strict versions this is only true statistically (e.g. Elitzur 
Additionally, Sy 1 galaxies should 



2006) 



& Shlosman 

have an unobscured point source (the unresolved accre- 
tion disk and inner rim of the torus) - but the relative 
strength of these two components in the mid-IR is model- 
dependent. 

To test this, we can compare our observations with 
other MIDI observations: In the C ircinus galaxy (Ltjv ~ 
4 x 10 36 W, i to ru S ~ 5 x 10 35 W), |Tristram et al | fl2007| 
found a warm (T ~ 330 K) disk with major axis FWHM 
~ 0.4 pc and a larger, similarly warm, component of ~ 2 
pc FWHM. 

In NGC 1068 (L uv ~ 3 x 10 37 W, L torus - 10 37 W), 



Raban et al.| ( p009| ) found a hot (T ~ 800 K) disk of 1.35 
x 0.45 pc and a warm component 3 x 4 pc in FWHM. 
They identified the disks with the densest parts of the 
torus of the unified model. 

For NGC 4151 (L uv - 1.5 x 10 36 W - variable (from 
NED), Ltorus ~ ^D 2 vF g - 6 x 10 35 W) we determined 
a torus size (FWHM) of ~ (2.0 ± 0.4) pc and a dust 
temperature of (285t^) K. 

When scaled to the accretion disk luminosity Luv, 
these values agree well with the torus sizes (r ~ Lyv 1 ^ 2 ) 
and torus luminosities (Ltorus ~ -^uv) of the two Seyfert 
2 galaxies and the temperatures are also very similar. 
Note, however, that the torus luminosity used here is 
calculated from the 12/xm flux density, not taking into 
account emission at longer wavelengths. 

4.2. Greybody models and the nature of the extended 

source 

In addition to the single-wavelength Gaussian models 
discussed above, we can also connect the various wave- 
lengths together by constructing greybody models where 
we assume a smooth powerlaw dependence of tempera- 
ture with radius: T(r) oc r~ a , and an essentially constant 
emissivity with radius and wavelength. With such mod- 
els we get acceptable fits (reduced x 2 ~ 1) to both the 
total and the correlated flux spectra with T(lpc) ~ 250 
K and 0.35 < a < 0.6. These values of a are consistent 
with dust receiving d irect radiation from a central source 
(e.g. |Barvainis|1987 1, i.e. an optically thin medium with 
optically thick clumps in it. 

These models can be extrapolated to shorter wave- 
lengths to check their consistency with the K band mea- 
surements. All plausible extrapolations of the MIDI data 
yi eld K band fluxes of < 1 mJy, much less th an observed 
by |Swain et al~j ( |2003[ ) and |Riffel et~aT| ( |2009|) (see below). 
On the basis of these greybody models, we therefore con- 
clude that the K band emission arises from structures 
which can probably not be extrapolated from the larger 
structure seen in the N band. 

From the greybody models one can further infer an 
emissivity e ~ 10 _1 - similar to what has been seen in 
NGC 1068 and Circinus. 

To conclude: The resolved nuclear mid-IR structure in 
NGC 4151 has a size, temperature and emissivity that is 
comparable to those in type 2 objects where the existence 
of clumpy tori is established. Apart from the similarities 
with type 2 tori and the temperature profile of the grey- 
body models, there is further evidence for dumpiness in 
the NGC 4151 torus from radio observations: IMundelll 



et al. (20031 measured HI absorption against the radio 
jet {PA ~ 77°) and found a structure of ~ 3 pc in size. 
From the velocities they further suggest that the gas is 
distributed in clumps. Warm dust possibly traces the HI 
gas in the mid-IR. 

It therefore seems reasonable to identify the warm dust 
structure resolved now in NGC 4151 with the clumpy 
tori seen in Sy 2 galaxies. Due to the limited observing 
geometry and the limited amount of observations we can- 
not reconstruct its apparent shape nor can we constrain 
a model with more than the two components discussed 
above. 

4.3. The point source and its relation to K band 
measurements 



K b and interferometry measurements by |Swain et al. 
(2003) and, with higher significance, more recent Keck 
observations (Pott et al. 2009, in prep.) revealed a 
marginally resolved source, compatible with hot dust at 
the sublimation radius of ~ 0.0 5 pc. Reve rberation mea- 



surements by |Minezaki et al. ( 2004 ) and |Koshida et al. 
(2009 1 find lag times At between the UV/optical contin- 
uum and the K band corresponding to a radius of ~ 0.04 
pc (variable) which they interpret as the sublimation ra- 
dius of dust; 

65 mJy at 



|Riffel et al.| ( 2009 1 measured a flux of 
2.2 /zm and from spectral fitting they found that, at 2.2 
/im, this flux is entirely dominated by a blackbody with 
a temperature of (1285 ± 50) K - again consistent with 
hot dust at the sublimation radius and the K band inter- 
ferometry, taking into account tha t L\jy is variable by 
at least a factor of 10 ( |Ulrich|2000[ ). To account for that 
variability when comparing our observations with these 
K band observations of different epochs, we looked at 
the X-Ray flux as a proxy for the UV-optical radiation 3 
and find that, at the date of our observations, the source 
was p r obabl y in a higher state than at the time of Riffcl 
et al. ( |2009[ )'s observations and emitted ~ 0.2 Jy in the 
K banal 

Since the flux density F v of a ~ 1285 K blackbody is 
roughly the same at 2.2/im as at 8.5 /zm, we can com- 
pare the flux in the K band (~ 0.2 Jy at the time of 
our measurement) with our point source flux at 8.5 /im 
(~ (0.119 ± 0.016) Jy). From this it seems likely that 
hot dust is contributing to our point source at 8.5 /im. 
The spectrum of the point source rises by more than a 
factor of two from 8.5 /zm to 12.5 /im, however, and thus 
requires emission from an additional small, "red" compo- 
nent. This could be core synchrotron emission, although 
such sources usually have flat spectra, or emission from 
a small, cool, optically thick central dust structure, pos- 
sibly shadowed from direct accretion disk radiation. 



5. CONCLUSIONS 

Using mid-IR interferometry, we have resolved a warm 
dusty structure in NGC 4151. Its FWHM size (2.0 ± 0.4 
pc - from comparing the data with a Gaussian model), 
temperature (285^50 K) and emissivity (~ 0.1) are in 
good agreement with the clumpy tori seen in type 2 

3 Data from the All-Sky-Monitor (ASM) on the Rossi 
X-Ray Timing Explorer (RXTE), available online at 
http://xte.mit.edu/asmlc/srcs/ngc4151.html The lag due to 
the light travel time ~ cAt ~ 48 days was taken into account. 
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AGNs and are thus consistent with the unified model 
of Active Galaxies. Excess emission around 10.5 /zm on 
the intermediate baseline indicates that silicate emission 
might be distributed out to ~ 1 pc in AGNs. 

Using simple models we compare our mid-IR fluxes 
with observations in the K band and find that the struc- 
ture we resolve is probably not the smooth continuation 
of the nuclear source detected in the K band 

Due to the limited number of measurements, no two- 
dimensional information could be gathered and questions 
about the unified model (such as: is the dust structure in 
Sy 1 galaxies thick and torus- like or rather disk-shaped?) 
remain unanswered. Since nuclear dust distributions are 
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